Zipping and collapse of diblock copolymers 
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Using exact enumeration methods and Monte Carlo simulations we study the phase diagram 
relative to the conformational transitions of a two dimensional diblock copolymer. The polymer is 
made of two homogeneous strands of monomers of different species which are joined to each other at 
one end. We find that depending on the values of the energy parameters in the model, there is either 
a first order collapse from a swollen to a compact phase of spiral type, or a continuous transition 
to an intermediate zipped phase followed by a first order collapse at lower temperatures. Critical 
exponents of the zipping transition are computed and their exact values are conjectured on the basis 
of a mapping onto percolation geometry, thanks to recent results on path-crossing probabilities. 
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I. INTRODUCTION 

Polymers in solution typically undergo a coil - globule 
transition from a high temperature (T) swollen phase to 
a low T phase where the polymer assumes compact con- 
formations. In the case of homopolymers, for which all 
the monomers are identical, this transition is by now well 
understood . It is known as 9-collapse and has been 
widely investigated in the past years using various meth- 
ods such as mean field approximations || , exact enumer- 
ations of interacting self-avoiding walks on lattices 
Monte Carlo ||, transfer matrix || and field theoreti- 
cal calculations In two dimensions the exponents of 
the 9-collapse have been related to the fractal properties 
of the percolation cluster and are believed to be known 
exactly f§. 

The study of the conformational properties and phase 
transitions of macromolecules with inhomogeneous or 
random sequences of monomers is an interesting frontier 
in nowadays polymer statistics These systems pose 
new theoretical and numerical challenges, compared to 
their more standard, homogeneous counterparts. Partic- 
ularly interesting is in general the possibility that the 
inhomogcncitics along the chain could lead to transitions 
and universality classes of scaling behavior, which are not 
realized for homopolymers |To|JTl[ ] . Moreover, the most 
complex versions of models of this class are also expected 
to be useful for the description of phenomena like protein 
folding [[ID, DNA denaturation Q and RNA secondary 
structure formation J3. Thus, an investigation of the 
universal properties of the simplest among these systems 
can offer an important gauge of the relevant model ingre- 
dients necessary in order to reproduce the basic confor- 
mational mechanisms in more sophisticated descriptions. 

One of the most elementary conformational transitions 
(not realized in homopolymers) one can try to describe 
in relatively simple terms is what we call here a zipping 
transition. By zipping we mean a process in which two 
strands composing the polymer come in contact in such a 
way as to form a bound double structure, which remains 



swollen and does not assume compact configurations. In 
order to induce a transition from unzipped to zipped 
state, the minimal inhomogeneity required implies a dis- 
tinction between the two strands: if the polymer is made 
of two blocks composed of monomers of different species 
(diblock copolymer) and there is a dominant attractive 
interaction acting between these different monomers, one 
would expect such a transition to be possible. Of course, 
the zipping occurring in biomolecules results in general 
from higher degrees of inhomogeneity than those of a 
simple diblock copolymer. 

From a physical point of view one can think to diblock 
copolymers with oppositely charged monomers in the two 
blocks; in the model discussed here the interactions are 
of short range, and this would correspond to the case 
of screened Coulomb forces. Another possibility is that 
the attractive interactions between monomers of the two 
blocks are established through a preferential formation of 
hydrogen bonds. The attractive interactions between the 
two blocks, besides zipping, tend also to produce collapse 
into a globular compact state, unless some contrasting ef- 
fect limits the capability of a given monomer to attract 
monomers of the other block. 

In a recent paper p5| , a model of diblock copoly- 
mer with some of the features discussed above has been 
studied in both two dimensions (2D) and 3D. In that 
model the two blocks were represented by two halves of 
a self avoiding walk (SAW) on hypercubic lattice, and at- 
tractive interactions were acting between nearest neigh- 
bor sites (monomers) visited by the two blocks. So, 
apart from the steric constraints, there was no inter- 
action mechanism possibly opposing the tendency of a 
given monomer to be surrounded by as many as possi- 
ble monomers of the other block. The transition of the 
diblock copolymer from a high T swollen to a low T com- 
pact phase, had analogies with both polymer adsorption 
on a wall, and 8-collapse, but turned out to belong to 
a universality class different from both p5[ . An intrigu- 
ing question remained open concerning the very nature 
of this transition: indeed, the possibility that a zipped, 
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swollen phase could exist for temperatures just below the 
transition, could not be excluded. If this were the case, 
the adsorption-like collapse found in Ref. jL{| would cor- 
respond to a zipping and a further transition to the com- 
pact globular phase should be expected to take place at 
a lower T. 

In the present article we extend the model of Ref. |l5| ] 
in 2D to include an interaction among alternating triplets 
of different monomers. Depending on its sign, this addi- 
tional interaction can enhance the tendency of the sys- 
tem either to form compact structures, or to take zipped 
conformations. We draw an accurate phase diagram for 
the system, in which a zipping transition line is well- 
identified and characterized. Our analysis seems to indi- 
cate that the adsorption-like collapse of Ref. belongs 
to the zipping universality class as well. Specifically, we 
find that, depending on the triplet interaction energy one 
has either a continuous swollen-zipped transition followed 
by a first order collapse into compact conformations, or 
a direct first order swollen-compact transition. Although 
we mainly focus on the zipping, it turns out that also the 
first order collapse has interesting features, since it shows 
remarkable analogies with that found in homopolymers 
with orientation dependent interactions, which attracted 
some attention recently 

We will argue that the exact exponents for the zipping 
transition in 2D can be found through an identification 
between the stochastic geometry of the blocks and that of 
a percolation cluster backbone p7| . A preliminary form 
of this argument was presented in Ref. ||15| . The relevant 
dimensions of the percolation cluster can be identified 
thanks to some recent results for path-crossing probabil- 
ities ]l8| . Our numerical estimates for the zipping expo- 
nents are in very good agreement with the conjectured 
values. The mapping onto percolative stochastic geome- 
try we discuss here is, to our knowledge, the first example 
of exact results derivation for a genuinely inhomogeneous 
polymer problem in 2D. The connection with percolation 
geometry shows also that the physics of zipping is closely 
connected with that of the ©-point transition, for which a 
representation in terms of percolation geometry has been 
established already long ago Q . 

Besides the prototypical importance that zipping ac- 
quires here also on the basis of our exact results, one 
should realize that many of the conformational transi- 
tions occurring in biomolecules show aspects which, to 
some extent, are reminiscent of zipping. This is certainly 
the case for DNA, in which, upon lowering temperature 
below the denaturation one, the conjugate bases form 
pairs, so that the molecule arranges itself in a double 
stranded helical structure jljj]. Double stranded, zipped 
structures appear also in the folding of /3-hairpin pep- 
tides JO)] . In the philosophy mentioned at the beginning 
of this section, one can hope that studies like the one 
presented here can teach something about how to model 
properly these more complicated systems. 

This paper is organized as follows: In Section || we 
present the model and the main features of the phase 



diagram. In Sections III and IV we discuss the numer- 



ical results obtained by exact enumerations and Monte 
Carlo simulations, respectively. In Section [y], using re- 
cent results for crossing probabilities of percolation paths 
in 2D, we conjecture the exact values of the exponents of 
the zipping transition. Section VI concludes the paper, 
with a summary of the results and a general discussion. 



II. THE MODEL AND PHASE DIAGRAM 



We model the 2D diblock copolymer by an interacting 
SAW on the square lattice. In the configuration w the 
SAW has |w| = N vertices (monomers), with N even, and 
consists of N/2 consecutive monomers of type A (wa) 
followed by N/2 monomers of type B (wg). A pair of 
vertices (A, B) form a contact if they are a unit lattice 
distance apart. The interaction between the two blocks 
wa and wb is taken into account by assigning an en- 
ergy e (e < 0) to each A — B contact. In addition we 
introduce a second energy parameter S, associated with 
contacts formed by a sequence A-B-A or B-A-B of neigh- 
boring monomers on a line. We refer to these sequences 
as to triple contacts. The Hamiltonian of the system in 
configuration w is given by: 



H(w) = N A b{w)s + N 3 (w) 6 



(1) 



where Nab(w) and N^(w) are the number of A — B and 
of triple contacts, respectively. 

For S — we recover the model introduced in Ref. [^5) ; 
in the present work we consider both signs of 5: a positive 
value of (5 must prevent the polymer from collapse into 
compact conformations and favor an intermediate zipped 
phase, while for a negative 5 the tendency to collapse is 
enhanced. 

Letting cn(Nab, A3) be the number of copolymer con- 
figurations with N edges, Nab contacts of type A — B 
and A3 triple contacts, we define the finite- A free energy 
per monomer 



F N (0,6) = N- 1 logZ N (p,5), 



where 



Z N {(3,5) = J2 c N (N AB ,N 3 ) e 

N A B;N 3 



/3(N AB s+N 3 S) 



(2) 



(3) 



is the partition function and /3 = 1/T |2Cf| . Throughout 
the rest of the paper we set e = — 1. By varying T and 
5 we have explored the phase diagram of the model (see 
Fig. |l|) on the basis of exact enumeration and Monte 
Carlo simulation results. 
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FIG. 1. Phase diagram in the 5-T plane. Line (1) is the 
continuous zipping transition line separating the high tem- 
perature swollen phase from the zipped phase. The lines (2) 
and (3) are of first order type. 

As expected, we find a zipped phase in the positive S 
region, while for negative 8 there is a direct transition 
from the swollen to the compact phase. The numerical 
results show that the line (1) separating the swollen and 
zipped phases is continuous, while (2) and (3) are first 
order. Unfortunately, the numerical methods at our dis- 
posal are not sufficiently accurate to determine the pre- 
cise location of the intersection point between the lines, 
nor the character of the phase transition in this point, 
which could be of special type. The location of this in- 
tersection point seems to fall at slightly negative values 
of 5. 

It is rather instructive to show some typical equilib- 
rium conformations of the copolymer for various values 
of T and S, in the three different phases (see Fig. |2|). The 
configurations are snapshots obtained by Monte Carlo 
simulations: (a) and (b) are conformations in the swollen 
high T phase, with (b) close to the zipping transition, (c) 
and (d) are instead zipped configurations, with (d) sam- 
pled in the vicinity of the transition to the compact state; 
notice that the pairing of the two strands in (d) follows an 
opposite orientation with respect to (c). Finally (e) and 
(f) are both compact, but of different nature: the latter 
occurs at 6 < where triple contacts are energetically fa- 
vored. Therefore the polymer assumes a spiral-like shape 
with straight segments turning around the center in order 
to maximize the number of triple contacts. In the case 
(e) 6 is positive: the configuration is still of spiral type, 
but in this case the arms of the spiral are oriented pref- 
erentially at 45° with respect to the axes of the square 
lattice, in order to avoid the formation of triple contacts. 

One can easily understand now why the lines (1) and 
(3) for large 6 run practically horizontal in the phase dia- 
gram of Fig. |l|. In the whole zipped phase triple contacts, 



which cost an energy <5, seldom occur (see Fig. ||(c)), 
therefore the zipping temperature should depend rather 
weakly on S. In addition the polymer can form compact 
conformations such as that shown in Fig. §(e), which 
also avoids this type of contacts. Hence, also the zipped- 
compact transition temperature should not depend on 6, 
when S is positive and large enough. Both lines (I) and 
(3) in Fig. j] should be asymptotically horizontal for large 
S. 




FIG. 2. Typical Monte Carlo equilibrium configurations for 
an N = 60 diblock copolymer: a-b) swollen chains; c-d) 
zipped chains; e-f) compact chains with S = 1.5 (e) and 
8 = -1.5 (f). 



III. EXACT ENUMERATIONS 

Exact enumerations of interacting SAW's are standard 
techniques for the study of the homopolymer O-collapse 
transition In the present calculation we generated 
all possible configurations for copolymers up to N = 30 
monomers, a length which is already sufficient to char- 
acterize rather well the critical behavior of the zipping 
transition. 

The occurrence of phase transitions in interacting poly- 
mer systems can be detected by studying the large N 
behavior of the canonical average squared radius of gy- 
ration 



Ri 



Z w exp[-PH(w)] R 2 (w) 



(4) 



where the sums extend to all A-steps configurations, w, 
of the copolymer, with radius R(w) relative to the center 
of mass. Indeed, in the proximity of a conformational 
transition temperature T c we expect that 



R g {N, T) ~ N Uc lZ [(T - T C )A*] 



(5) 
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where v c and (f> are the exponents characterizing the tran- 
sition and 1Z is a scaling function that is assumed to 
approach a positive constant if its argument approaches 
zero. 

Another important quantity is the specific heat 
C{N,T) = ±d{H)/dT, which for N large and T close 
to T c is expected to obey the scaling: 



R'(N,T = T C ) 



C(N, T) ~ N^^C [{T - T C )A*] 



(6) 



where C is again a suitable scaling function. 

Figures || and [| show dR 2 /d(3 and C as functions of f3 
for S — (a) and 6 = — 1.5 (b). As the radius of gyra- 
tion drops at a transition, its derivative dR 2 /df], shows 
a peak in correspondence to the transition point. 

For S = —1.5 both quantities have a single isolated 
peak, indicating that there is a single transition from 
a swollen to a compact phase. For S = 0, instead, the 
derivative of the radius of gyration has two distinct peaks 
(Fig. 0(a)), while the picture emerging from the specific 
heat (Fig. ||(b)) is somewhat more confusing, since the 
iV-dependence of the peak positions and heights is rather 
irregular and their extrapolation to N — > oo becomes im- 
possible. Therefore we focus on the peaks of dR 2 /d(3. Ex- 
trapolating their positions at 5 = we finds the following 
two estimates of critical temperatures T c = 1.40(15) and 
Tic = 1-1(20), which overlap somewhat within error bars. 
For this reason it is difficult to discern between two sepa- 
rate, but close, transitions, and a single one. However, as 
S is increased, the two sets of peaks get clearly separated 
and extrapolations yield two distinct transition temper- 
atures. We focus here on the characterization of the high 
T transition from the swollen to the zipped phase. The 
scaling form of Eq. (j^) implies that T C (N), the temper- 
ature where dR 2 /d/3 has a maximum, scales for large N 
as T C (N) ~ T c ~ xqN^^, with xo a suitable constant. We 
calculated both the radius of gyration and specific heat 
at T C (N); from (§) and (g) one has: 



and 



R g (N,T = T c (N))~N"° TZ(x ) 



C(N,T = T C (N)) ~ N 2 ^ 1 C{x ) 



(7) 



(8) 



For the calculation of the critical exponents we formed 
first the finite N approximants, for instance 



HR g (N + 2)/R g (N)) 
ln((JV + 2)/N) 



(9) 



(here R g (N) is a shorthand notation for R g (N,T c (N)j) 
and then extrapolated v c (N) to N — > oo. The same pro- 
cedure was followed for <j>. 

The extrapolated values are reported in Table |, to- 
gether with the exponent u' c , which is that associated to 
the radius of the half-chain, or single block, which at the 
critical temperature should scale as: 



(10) 



TABLE I. Extrapolated values of v c and <j> from the ex- 
act enumeration data relative to diblock copolymers up to 
length TV = 30. The exponent v' c is obtained from the scaling 
behavior of the radius of one of the blocks. 



s 


Vc 


<$> 




0.0 


0.72(1) 


0.60(5) 


0.74(1) 


0.5 


0.73(1) 


0.58(3) 


0.750(5) 


1.0 


0.73(1) 


0.57(3) 


0.750(5) 


3.0 


0.74(1) 


0.56(3) 


0.750(5) 




FIG. 3. Derivatives of the squared radius of gyration with 
respect to the inverse temperature for 5 = (a) and S = —1.5 
(b) and for TV = 12, 14, . . . 30. The double peak structure for 
S — suggests the sequence of two transitions, swollen-zipped 
and zipped-collapsed. 

The values of the exponents <fi and v c vary slightly 
along the line (1) when S is increased, while v' c is rather 
stable. We believe that the variation of v c and <f> is a spu- 
rious effect due to the vicinity of an additional transition 
in the neighborhood of 5 — 0. It is much more plausible 
that the exponents are constant along the line (1); the 
most reliable estimates for v c and <p should be those for 
large S, where the lines (1) and (3) are clearly separated. 
The values of v c and v' c are consistent (the former only 
at large 8) with the scaling behavior of a SAW, namely, 
R g ~ A 3 / 4 . The value of <fi is instead consistent with 
4> = 9/16 = 0.5625, which was conjectured in Ref. (T^j 
for the transition at S — and will be derived in detail 
m Section 0. 

As for the transition lines (2) and (3), the exact enu- 
meration analysis is not at all conclusive since the scaling 
behavior of the peaks with the chain length N is rather 



4 



irregular and precise extrapolations turn out to be impos- 
sible. This issue will be clarified with the use of Monte 
Carlo simulations, which allow to achieve much larger 
copolymer lengths. 




FIG. 4. Solid lines: Specific heat from exact enumerations 
for N = 12, 14, ... 30 for S = (a) and S = -1.5 (b). Circles: 
Monte Carlo results for N = 30. 



and B paired together, but still not compact (see Fig. [T] 
(b,c,d)). To increase the mobility of the Markov chain in 
this region we added a set of bilocal moves, such as end- 
end reptation and kink-end (and end-kink) moves [B7) . 
These moves are particularly effective for dense chains, 
and even more effective for the zipped chains, where typi- 
cally one side of each half-chain is free and can hold a new 
kink. The resulting algorithm is a little heavier, but en- 
ables the reciprocal sliding of the half-chains and a more 
efficient exploration of the configuration space. One may 
then run in parallel a number m (typically 20 — 40) of 
these Markov chains at different temperatures. The sam- 
pling at low T is then considerably enriched by swapping 
configurations between Markov chains contiguous in T. 
The whole process is itself a (composite) Markov chain, 
obeys detailed balance and is ergodic pi] . 

Monte Carlo simulations were performed for three dis- 
tinct values of 5: 6 = 1.5, 5 = —1.5 and 5 — 0. As a test 
of the performance of the multiple Markov chain algo- 
rithm we compared the Monte Carlo results with those 
obtained from the exact enumeration for chains up to 
N = 30 monomers. In all cases analized the agreement 
turned out to be extremely good (see, for instance, Fig. 
|(a)). 



A. 5 = 1.5 



IV. MONTE CARLO SIMULATIONS 

In order to sharpen and extend the results of exact enu- 
merations, in particular concerning the properties of the 
zipped and collapsed phases, we performed Monte Carlo 
simulations for various N and 5 and for a wide range of 
temperatures. Since the simulations considered involve 
sampling at points which include low values of T, a stan- 
dard Markov chain Monte Carlo approach is unlikely to 
be successful, being difficult to construct a Markov chain 
sufficiently " mobile" at low T where the interaction ener- 
gies become relevant. Instead we use a multiple Markov 
chain technique by which one samples simultaneously at 
various values of T, including T = oo where converge is 
rapid. Most recently this method has been used success- 
fully to investigate collapse transitions in homopolymers 
[21 1, hetheropolymers pdj ], and adsorption in 0-solvent 

First one defines a Metropolis based Markov chain for 
a temperature T. This procedure makes use of a hy- 
brid algorithm based on pivot |24| as well as on local 
moves p5| . Pivot moves are of global type and oper- 
ate well in the swollen regime, whereas local moves turn 
out to be essential in speeding up Monte Carlo conver- 
gence at low temperatures [^l). In our calculations, each 
Monte Carlo step consists of 0(1) pivot moves and O(N) 
local moves. In this model, however, we have to deal 
also with a zipped phase where the most probable con- 
figurations are characterized by having the two blocks A 



In the case 5 = 1.5 we considered diblocks of lengths 
up to N = 400 and sampled at a set of m ~ 40 temper- 
atures typically in the range T G [0.5, oo]. In Fig. || we 
plot the specific heat as a function of (3 for different N 
values. Clearly each curve displays a double peak struc- 
ture indicating two subsequent transitions. We can rule 
out the possibility that such double peaked structure is 
a finite size effect by noting that the peaks sharpen and 
grow with N. Let us focus first on the set of peaks at 
higher temperatures, i.e. on the transition from a swollen 
to a zipped phase. The corresponding T c and 4> could be 
deduced from the N dependence of the height, h(N), and 
position, T C (N), of the peak maxima. Indeed, from the 
scaling behavior (11) we expect, as N increases, 



h(N) ~ N 2 ^ 1 and T C (N) 



N -<t> 



(11) 



Since a linear least squares fit of log ft. vs. log TV gives 
a very large x 2 statistical error, we fit the data with 
a function AN 2 ^~ 1 {1 + B/N) where a scaling correc- 
tion 1/N is included. The least squares fit in this case 
gives 4> = 0.57 ± 0.02, in agreement with the value 
4> = 9/16 = 0.5625 conjectured in Ref. and also 
with the estimates obtained by exact enumeration. This 
procedure yields results consistent with a direct extrapo- 
lation of effective finite N exponents. We have also tried 
to fit the data by using the more general scaling cor- 
rection 1/N , but the best fit is obtained for A = 1. 
The estimated value of (p allowed us to extrapolate T c 
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by plotting T C {N) vs. I/AT*. This gives T c = 1.51(4) 
{fi c = 0.66(2)). 
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FIG. 5. S — 1.5: Plot of the specific heat vs. /3, for 
chains of various lengths. For N = 300 and N = 400 the 
low temperature peak was not reached because of the con- 
siderable autocorrelation time caused by low mobility of long 
collapsed chains. Inset: collapsed specific heat C//V 2 ^ -1 vs. 
(/3 - AOJV*, with <j> = 9/16 = 0.5625 and (3 C = 0.66. 

The inset of Fig. | shows a plot of C/N 2 ^- 1 vs. 
(fi — P C )N^, where we have used <fi = 9/16 and the es- 
timate p c = 0.66. As expected, the high temperature 
peaks collapse onto a single curve quite nicely. On the 
contrary, for the set of peaks at lower temperatures, the 
same rescaling procedure turns out to be inappropriate. 
In particular, by using the rescaled variables adequate 
for the former set of peaks, the positions of the latter set 
tend to move away from zero while their heights still in- 
crease with N: this is only consistent with a scenario in 
which a new transition, at a lower T = T2 C , exists. This 
transition should be also characterized by a crossover ex- 
ponent greater than <fr = 9/16. We have tried to verify 
this by looking for two new values /?2 C and 4*2 that al- 
low a reasonable fit of the scaling behavior of the second 
set of peaks. In this way we obtained the rough esti- 
mates p2c ~ 1.0 and cf>2 ~ 0.7 > 9/16. Unfortunately 
the sampling at low temperatures is not good enough to 
make such estimates sufficiently sharp. Moreover, in the 
zipped phase the effective size of an N monomers system 
drops to N/2, making the finite size corrections to scaling 
more pronounced. 

The different nature of the two transitions can be bet- 
ter detected from the behavior of P(E,N), the proba- 
bility distribution of the energy E, for a chain of length 
N. Figure [| shows a plot of P(E, N) as a function of 
E/N for N = 200. At sufficiently high temperatures this 
quantity has a maximum in E = and decreases rapidly 
with E/N (Fig. ^(a)). As the temperature is lowered 
the maximum shifts continuously to larger values of E/N 



(b,c). For lower temperatures P develops a double peak 
structure (d,e,f). In the case (e) the peaks have equal 
height, while at temperatures below or above it one of 
the two peaks dominates over the other. This behavior, 
which persists and becomes more pronounced upon in- 
creasing N , is an indication of phase coexistence; hence 
the transition at lower T is of first order type. In terms 
of specific heat this would mean C(T) ~ N as T = T2 C , 
i.e. (f>2 — 1. By extrapolating /?2c(A r ) vs. 1/AT we find 
foe = 1.2 ±0.2. 
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FIG. 6. Plot of P(E, N) with TV = 200 and for various 
temperatures: (a) (3 = 0.3, (b) (3 = 0.75, (c) (3 = 1.4, (d) 
(3 = 1.58, (e) (3 = 1.66 and (f) (3 = 1.8. 

Another way to characterize the different phases of the 
model consists in looking at the scaling behavior of met- 
ric quantities such as R 2 g defined in (Ji|) and the mean 
squared end-to-end distance R 2 (AT) = ((rjv — r ) 2 }, where 
ro and rjv are the two end-monomers of the copolymer. 
For large A" we expect 

R 2 e (N) ~ Pe N 2 » (12) 
R 2 g (N) ~ p g N 2 » (13) 

and an interesting quantity to be computed is the ratio 
Pe/Pg, which is expected to be universal |23). For non in- 
teracting SAW's on a square lattice, exact enumerations 
and Monte Carlo simulations give p e /p g ~ 7.13 (see ref- 
erences in p3|). Figure ^| shows the ratio R 2 /R 2 as a 
function of (3 for several AT. Note that in the range of 
< (3 < 0.66 the curves tend to assume a constant value 
Pel p g — 7.15(5) in agreement with the value expected 
for non interacting SAW's. In the proximity of the zip- 
ping transition the curves start to bend downwards and 
at 0c ~ 0.66 they cross each other almost in a unique 
point (see inset). At the crossing point our estimate of 
the universal amplitude is p e / p g = 6.35 ± 0.20 which 
is definitively different from the amplitudes ratio of the 
SAW universality class. The zipping transition that can- 
not be distinguished from the swollen phase in terms of 
the v exponent, is however characterized by a different 







value of the universal ratio Pel Pg 22 




FIG. 7. R e /R g vs. fi for 5 = 1.5 and for various chain 
lengths. The N = 30 curve is calculated from exact enu- 
merations, while the others are obtained from Monte Carlo 
simulations. Inset: blow up of the crossing region. 

If the temperature is further lowered the R 2 /R 2 g curves 
reach a minimum value that decreases as N increases. 
This is an indication that the end-to-end distance in the 
zipped phase does not scale anymore like the radius of 
gyration, as assumed in Eqs. (|l2|,|l3|). For sufficiently 
low temperatures, p e /p g starts to grow back indicating 
that the compact phase is characterized by end-to-end 
distance and mean radius of gyration that scale in the 
same way with N . As the typical low T configurations 
are of spiral type (see Fig. ||(e,f)) with end-points at 
opposite sides of the spiral, it is natural to expect that 
R e -R^N 1 ' 2 . 
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For 8 < triple contacts are favored and we expect 
(as the exact enumerations already indicate) a single 
transition from swollen directly to compact phase. To 
investigate the nature of such transition we have per- 
formed runs with S = —1.5, for several values of N, sam- 
pling at m ~ 30 different temperatures in the interval 
T G [1.3, oo]. As in the case 5 — 1.5 we have examined 
the probability of finding the copolymer in a configura- 
tion with energy E, as a function of the temperature. A 
plot of P for N = 200 and three different temperatures 
is shown in Fig. [|. Close to the transition temperature 
P has two maxima (see Fig. |§|(b)) one at E = and 
the other at E/N ~ 0.7. This is a clear indication of a 
first order transition. The evidence of such a behavior is 
stronger than in the case S = 1.5 since here the coexis- 
tence is between two phases (the swollen and the com- 



pact) with a rather large difference in energy, therefore 
the double peak structure of P can be noticed already 
for small N. 

From the analysis of the specific heat peaks we find 
that they become sharper as N increases and their height 
appears to grow with a power of N slightly exceeding one 
(the physical upper limit). At the same time, in a plot of 
energy vs. temperature we see curves that seem to ap- 
proach step functions. These data support the idea that 
the corresponding transition should be of first order. 
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FIG. 8. P(E,N) vs. E/N for S = -1.5 and N = 200 for 
three different temperatures (a) (3 — 0, (b) (3 = 0.545 and (c) 
(3 = 0.6. Inset: Plots of P(E,N) vs. E/N near phase coex- 
istence for TV = 60, = 0.66 (solid line), N = 100, (3 = 0.56 
(dotted line) and N = 140, (3 = 0.545 (dashed line). 



C. 8 = 

The more delicate region to be explored in the phase 
diagram is the neighborhood of 5 = 0, where three transi- 
tion lines meet each others. We have chosen in particular 
the case 6 — 0, since it was considered in Ref. ]l5[ |. 

In Fig. H we plot the specific heat as a function of (3 for 
several N values. For the smallest chains (N = 60,80), 
one observes a peak in the specific heat with a shoulder 
at smaller (3. When the copolymer length is increased 
the shoulder becomes hardly noticeable. From the spe- 
cific heat plot one cannot rule out the possibility that the 
shoulder eventually vanishes leaving out a single transi- 
tion from a swollen to a compact phase. The other possi- 
bility is that there are two separate transitions, but very 
close in temperature. 

The presence of two distinct transitions is suggested 
by a plot of the temperature derivative of the total ra- 
dius of gyration, shown for N = 100 and N = 200 as 



thick lines in Fig 
two peaks, which a 



I0[ In this case one clearly detects 
though coming closer to each other 
as N increases are still noticeable and sharp for N rather 
large. The thin lines in Fig. 00 are the temperature 



7 



derivatives of the radius of gyration of a single block, 
which show only one peak in correspondence to the low 
temperature peak of the derivative of the total radius of 
gyration. This behavior is consistent with the following 
picture: coming from the swollen phase (small (3) one has 
first a zipping transition characterized by a drop of the 
total radius of gyration, while the radius of gyration of 
a single block still behaves as a SAW and is not sensible 
to the zipping transition. However, at lower tempera- 
tures, in correspondence to the transition from zipped to 
compact phase both quantities drop and their derivatives 
show a peak. 

1.5 | 1 




P 

FIG. 9. 8 — 0: Plot of the specific heat vs (3, for various 
chain lengths. 
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FIG. 10. 8 = 0. Thin lines: plot of the derivative of the 
half-chain squared gyration radius, as function of (5. Thick 
lines: two examples of derivative of the total-chain squared 
gyration radius (N=120, 200). 

Another quantity which we investigated is the univer- 
sal amplitude ratio between the end-to-end distance an 



radius of gyration squared, which is plotted in Fig. |llj. 
Here, as for the 6 = 1.5 case, this universal quantity takes 
the SAW value ~ 7.13 at high T and drops in correspon- 
dence of the transition. The fact that we find intersec- 
tions with R^/Rg w 0.635 (the same value as for 6 = 1.5) 
strongly suggests the presence of a zipping transition with 
the same universal properties as that at 6 = 1.5. Unlike 
from Fig. [j] here R\j ' R 2 g drops and increases again in a 
narrow range of (3 values indicating the the zipped phase 
is restricted to a small temperature interval. 

In summary, although the numerical evidence is not 
fully conclusive, our data seem to favor the existence 
of two separate transitions for 6 — 0. As in the case 
S = 1.5, it is natural to expect that the low T one (zipped- 
collapsed) is of first order type. 
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FIG. 11. 5 = 0: Square end-to-end distance over square 
gyration radius as a function of (3. 



V. PERCOLATION PATHS AND EXACT 
EXPONENTS OF THE ZIPPING TRANSITION 

In this section we present a conjecture on the rela- 
tion between the statistics of some percolation paths at 
threshold and the diblock copolymer zipping transition. 
This conjecture leads to predict exact values for the ex- 
ponents. A preliminary, less precise version of the argu- 
ments below was given in Ref. fll5| |. 

It is well-known that, in 2D, the statistics of a ring 
polymer at the 0-transition is identical to that of the 
external perimeter, or hull, of a percolation cluster. 
Through this identification the exact exponents of the 
0-transition, vq — 4/7 and (b® — 3/7, were derived [||. 
Here, we show how similar arguments can be invoked 
for the zipping transition. The differences are mainly 
associated to the fact that the relevant percolative set 
appropriate for the zipping is not the hull, as for the 
homopolymer 0-point. 

Like in the 0-point case, here it is convenient to con- 
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sider site percolation on a triangular lattice. For this 
problem the relevant percolation contours, like the hull 
of a cluster, are in fact strictly self-avoiding paths on 
the dual, hexagonal lattice. Thus, also the equivalent 
diblock copolymer problem realized here by percolation 
paths will be on hexagonal, rather than square lattice. 
On the basis of universality we expect our results to ex- 
tend also to the square lattice case. 




FIG. 12. Percolation cluster of connected occupied 
hexagons (dashed). Each hexagon is centered on a site of 
the dual, triangular lattice. 

Let us consider a percolation cluster as sketched in Fig. 
|l2| . Its external perimeter is a self-avoiding ring. The en- 
semble of all possible conformations of an external hull on 
the lattice can be regarded as a problem of ring polymer 
(grand canonical) statistics, as discussed in Refs. One 
further realizes that this effective ring polymer problem is 
characterized by attractive interactions. These originate 
from the fact that, at threshold (p — p c = 1/2), multiple 
visitations |2^| by the hull of the same occupied, or va- 
cant, hexagon, give higher probability to the realization 
of a ring configuration. Indeed, when the contour pro- 
ceeds essentially straight, to each new step corresponds 
a new hexagon whose state (occupied or vacant) has to 
be determined. So, each step implies a factor p c = 1/2 in 
the probability weight of the whole configuration. When 
the contour folds on itself and revisits, after some steps, 
the perimeter of the same hexagon, the factor 1/2 does 
not apply, resulting in higher global probability. This is 
equivalent to an attractive interaction favoring the mul- 
tiple visitations of the same hexagon. 

It is convenient here to summarize some very recent 
exact results concerning the fractal dimensions of vari- 
ous percolative sets. Following Rcf. [jl8| we consider an 
annular region of the hexagonal lattice delimited by an 
inner circle of small radius r, and an external one, of ra- 
dius R 3> r. Two types of paths connecting the two cir- 
cles are also considered. These paths are formed by con- 
nected and self-avoiding sequences of either occupied, or 
empty hexagons. The so-called path-crossing probability, 
namely the probability that I non-overlapping paths con- 
nect inner to outer circles, was found to behave asymp- 
totically as: 



Pi(r,R) « (r/i?)* 1 (14) 

where: 

I 2 - 1 

*« = — (15) 

The formula is valid if there is at least a path of each type, 
and the probability depends only on the total number of 
these paths, not on their type ]lq|. 




(a) (b) / 

FIG. 13. Path crossing configurations for (a) a dashed and 
a solid line and (b) two dashed and two solid lines. The proba- 
bilities of the configurations yield the fractal dimensions of the 
external perimeter of the hull (a) and of the cutting hexagons 
of the backbone (b). The dashed region indicates the perco- 
lating cluster of occupied hexagons, while the double dashed 
region of (b) shows a dangling end, a part of the cluster which 
does not belong to the backbone. 

Figure |l^ shows examples of crossing paths; without 
resolving the underlying lattice structure, we draw them 
as solid lines if they connect filled hexagons, while dashed 
lines are used for paths connecting empty hexagons. As 
a first example (see Fig. |l3|(a)) we consider the crossing 
probability for a continuous and a dashed path, which 
according to Eqs. (|T^ , [l5| ) decays as Pi = i ~ (r/i?) 1 / 4 . 
One recognizes immediately that the set of points for 
which two of such self-avoiding paths can be drawn are 
those of the external perimeter or hull of the percola- 
tive cluster (see Fig. |l^(a)). This identification al- 
lows to derive the fractal dimension of the hull. Since 
the area enclosed by the annulus is proportional to i? 2 , 
the perimeter of the hull enclosed in the annulus must 
scale as L eh ~ R 2 Pi =2 ~ R 2 ' X2 = R 7/4 - Identify- 
ing the external hull as a polymer ring at the O-point 
one then derives that the latter has a fractal dimension 
D e = A=2 = 2 - x l=2 = 7/4. 

In order to make contact with the diblock copolymer 
zipping, let us now imagine to identify two points, 1 and 
2, dividing the cluster hull in two equally long parts 
(see Fig. [14]). By fixing these two points on the clus- 
ter perimeter, one automatically defines a backbone as a 
subset of the whole cluster. The backbone is the union 
of all connected paths of occupied hexagons, which are 
strictly self-avoiding (i.e. in each path a given hexagon 
appears at most once) and join the points 1 and 2. In 
force of the definition, the backbone does not include 
the so called "dangling ends", i.e. those branches of the 
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cluster connected to the main body by narrow bridges 
(i.e. by regions in which only one occupied hexagon is 
available, making impossible for a self-avoiding path of 
hexagons to penetrate and exit at the same time). An 
example of dangling end is also schematically shown as 
double dashed area in Fig. |l3|(b). 




FIG. 14. (a) Schematic representation of a percolation clus- 
ter with dangling ends (dashed areas). When these are elimi- 
nated one remains with the cluster backbone (b) . The dashed 
segments cut the cluster in correspondence to the so-called 
"red" hexagons. 

The two points we hx on the contour clearly divide into 
two sides the perimeter of the backbone. Even if in this 
case it is not possible to give a simple expression of the 
effective interactions determining the shape of the two 
backbone sides, we expect them to be essentially local, 
like in the case of the hull, and to act differently ac- 
cording to whether they involve close encounters of the 
same side, or between different sides. This is consistent 
with the idea that the two sides of the backbone perime- 
ter could represent the statistics of a ring version of the 
diblock copolymer at the transition, the two parts corre- 
sponding, respectively, to blocks A and B. 

For the calculation of the fractal dimension of the ex- 
ternal perimeter of the backbone one can use Eqs. ( |il||l5| ) 
taking two continuous paths and a dashed one. This 
configuration clearly identifies the perimeter of the back- 
bone. Indeed, the two continuous paths guarantee that 
occupied hexagons inside the interior circle belong to a 
whole path connecting two infinitely distant points. At 
the same time, a dashed path implies that the vacant 
hexagons facing the occupied ones belong to the exterior 
of the cluster, and thus, also of its backbone. 

Therefore we now take I = 3 for the exponents defined 
in Eq. ( 15 ) . In this case we find that the external perime- 
ter of the backbone scales as ~ R 2 ~ Xi = i? 4 / 3 , which 
implies a fractal dimension D% = 4/3. This dimension is 
consistent with that found for the diblock copolymer at 
the zipping transition. Furthermore, it is natural to asso- 
ciate the switching on of effective attractive interactions 
between the two backbone sides to the existence of nar- 



row bottlenecks in the backbone itself (corresponding to 
only one hexagon). These are the so-called cutting or 
"red" hexagons of the backbone jl7|, which are visited 
by the two blocks simultaneously. In oder to determine 
their fractal dimension one has to consider a percolative 
configuration with two continuous and two dashed paths 
joining the circles, as sketched in Fig. |l3|(b). These iden- 
tify a dimension D4 — 2 — a; 4 = 3/4. Thus, for a back- 
bone with external perimeter equal to N and an average 
of Nab contacts we find N ~ R° 3 and N AB ~ R D K 
Consequently, the average number of contacts between 
the backbone sides grows like N AB ~ N Di /° 3 = N 9 / 16 . 
By identifying the external perimeter of the backbone 
with the ring diblock copolymer at its transition, one 
eventually finds v c = 3/4 and <j> — 9/16. Now the numer- 
ical determinations of v c and <j) at the zipping transition 
are remarkably consistent with these values, making the 
conjecture extremely plausible |3(l]. 



VI. CONCLUSIONS 

In this paper we studied the phase diagram for the col- 
lapse transition of a diblock copolymer with attractive 
interactions between monomers of different species and a 
triple contact interaction S, which according to its sign 
may either favor, or unfavor, compactification. In the 
region of negative S we find a first order transition from 
a swollen to a compact, spiral phase, while in the posi- 
tive 5 region there is a sequence of a continuous zipping 
transition and a collapse of first order type to compact 
conformations at a lower temperature. 

Our exact enumerations and Monte Carlo simulations 
yield numerical estimates of the critical exponents u c and 
<j> of the zipping transition, which are consistent with 
those we could conjecture using recent results for the 
fractal dimensions of the percolation cluster backbone, 
from which we expect v c = 3/4 and <j) — 9/16. The nu- 
merically determined exponents, therefore, are strongly 
supporting the hypothesis that the transition admits a 
description in terms of percolative stochastic geometry: 
the two blocks of the copolymer have the same fractal ge- 
ometry as the two sides of a cluster backbone, and their 
contacts correspond to the cutting hexagons, or links of 
the same backbone. This is, to our knowledge, the sec- 
ond example of a percolative representation for a poly- 
mer conformational transition in 2D, besides that of the 
0-point. The common percolative roots of these transi- 
tions suggests the possibility of a deep link between them, 
which ought to be elucidated by further studies. 

The results obtained for the various transitions appear- 
ing in the phase diagram help in clarifying the nature 
of the adsorption-like collapse occurring at 8 = and 
first detected in Ref . . In spite of the fact that most 
tests are not able to put into clear evidence the existence 
of two successive transitions, the only multicritical be- 
havior which can be characterized coming from the high 
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temperature region seems to definitely belong to the uni- 
versality class of the continuous zipping transition iden- 
tified for positive values of 8. Besides the compatibility 
of the exponent estimates, a very strong support to such 
conclusion comes from our determination of the univer- 
sal amplitude ratio between squared end-to-end distance 
and radius of gyration of the polymer. 

Another interesting aspect of the phase diagram cal- 
culated in this paper are the first order swollen-collapsed 
and zipped-collapsed transitions found, respectively, for 
negative and positive values of 8. In particular the lat- 
ter resembles the transition from swollen to spiral state 
found in oriented polymers |jl6|| , i.e. chains to which an 
overall orientation is assigned and where different ener- 
gies are associated to contacts between parallel or an- 
tiparallel segments of the chain. In fact the analogy be- 
tween the diblock copolymer in a zipped state and an 
oriented polymer is very appropriate: in the zipped di- 
block parallel contacts are of AB type, while antiparallcl 
ones are contacts between monomers of equal type. Dif- 
ferent energies are clearly associated with the two types 
of contacts. 

It is worthwhile to recall that simple polymer mod- 
els with some sort of zipping transition attracted already 
some attention in the recent literature pUB^ , mainly 
because of the relevance that such transition can have 
for biopolymers. Imbert et al. (3^] considered a diblock 
formed by two strands of oppositely charged monomers 
interacting with each other through long range Coulomb 
forces and found evidence of the existence of a zipping 
transition followed by a collapse at lower temperatures. 
Causo et al. |j2[ considered a simple model for the DNA 
denaturation transition, in which only the monomers 
which are at equal distances along the sequence from the 
center of the chain interact. They found evidence of a 
first order transition, from a swollen to a zipped phase. 
By its construction their model has no other transitions 
to compact state. In their case the first order zipping 
seems to be due to the selective interactions of monomers 
along the chain. Also in our model, if we turn on interac- 
tions only between AB monomers at equal distance from 
the center we find evidence of a first order zipping tran- 
sition. 

Finally, we point out that there are several possible ex- 
tensions of this work. First of all, it would be interesting 
to generalize the model to three dimensions and to inves- 
tigate the properties of the zipping transition in that case 
p5[ . Another open issue is the effect of disorder on the 
interaction between monomers for the zipping transition, 
which would allow to understand the behavior of models 
of polymers more relevant for applications to chemistry 
or biology than a simple diblock. 

We thank F. Seno for discussions and collaboration 
in the early stages of this work. Financial support by 
MURST through COFIN 1999 and INFM through PAIS 
1999 is gratefully acknowledged. 
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